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摘要 : 利用 组 成 星座 的 小 卫星 ， 分 别 携带 分 离 的 子 望远镜 和 合 光 成 像 望远镜 ， 构 成 
Fizeau 型 光学 综合 孔径 干涉 系统 ， 实 现 高 分 辩 率 的 面 源 目 标 成 像 是 当前 的 研究 热点 之 一 。 这 
种 光 干 涉 成 像 系 统 ， 由 于 稀疏 度 较 大 ，UV 覆盖 不 全 ， 即 空间 频率 采样 不 连续 ， 表 现 为 系统 
光学 传递 函数 有 零 值 存在 。 要 克服 UV 覆盖 不 全 的 影响 ， 获 得 等 效 的 大 孔径 望远镜 成 像 效 
I 果 ， 需 要 改变 子 孔 径 的 空间 排 布 ， 获 得 不 同 基线 条 件 下 的 图 像 ， 进 行 空间 频率 信息 的 提取 和 
Y 合成 ， 最 后 采用 着 滤 波 的 方法 达到 提高 图 像 质量 的 目的 。 在 分 析 单 子 孔 径 传 递 函数 与 系统 传 
递 函 数 关系 的 基础 上 ， 优 化 子 孔径 的 排 布 方式 ， 采 用 不 同 的 频 域 滤波 器 将 不 同 基 线 获 得 的 图 
像 中 信 嗓 比 高 的 频率 区 域 提 取出 来 进行 合成 ， 再 变换 到 空域 并 进行 着 滤波 处 理 ， 得 到 改善 的 
合成 图 像 。 仿真 结 果 显 示 ， 当 得 到 的 干涉 图 信 品 比较 低 时 ， 该 方法 可 以 有 效 地 提高 合成 图 像 
T 的 质量 。 
T 关键 词 : 稀疏 ;孔径 图 像 ;， 复原 频 域 滤波 ; 图 像 合成 
i 中 图 分 类 号 : TN911.73 文献 标识 码 : A 文章 编号 : 1672-7673(2016)03-0351-07 


由 于 孔径 的 衍射 效应 ， 传 统 望 远 镜 的 成 像 分 辨 率 受 限 于 人 瞳 口 径 。 受 加 工 工 艺 和 成 本 等 因素 的 制 
约 ， 光 学 系统 的 口径 不 可 能 无 限 增 大 "1 。 为 了 解决 这 个 问题 ,研究 人 员 提 出 了 光 干 涉 的 方法 。 随 着 
技术 的 发 展 ，Fizeau 型 光 干 涉 技 术 因 具有 一 定 的 成 像 视 场 和 对 面 源 短 时 成 像 的 特点 ， 受 到 越 来 越 多 的 
© 重视 ， 成 为 当前 的 研究 热点 之 一 。 早 在 1996 年 ， 美 国 空军 武器 实验 室 就 启动 了 用 于 建立 对 地 观测 的 
一 可 展开 式 望 远 镜 阵列 UltraLITE H H 
一 结合 近年 来 蓬勃 发 展 的 小 卫星 星座 技术 ， 人 们 提出 了 基于 小 卫星 的 Fizeau 型 光 干 涉 成 像 系统 。 
基本 方法 是 ; 在 每 个 小 卫星 上 各 自 携 带 一 台 小 口径 望远镜 ， 在 中 心 的 主 卫 星 上 携带 合 光 成 像 望远镜 ， 
子 卫星 和 主 卫 星 在 空间 上 按 一 定 的 位 置 排 布 ， 来 自 各 子 望远镜 的 光束 在 合 光 成 像 望 远 镜 的 焦 面 上 相干 
AM, 干涉 成 像 。 这 种 形式 的 光 干 涉 成 像 系统 ， 受 小 卫星 避 撞 设计 的 需要 ， 以 及 小 卫星 有 效 载荷 质量 
和 体积 的 限制 ， 子 望远镜 的 口径 和 要 实现 的 衍射 极限 分 状 率 所 对 应 的 口径 相 比 ， 必 然 差 别 很 大 ， 属 于 
稀 玻 度 很 大 的 光学 综合 孔径 系统 。 此 时 ， 由 于 UV 覆盖 不 全 ， 系 统 不 可 避免 地 对 相应 的 空间 频率 在 采 
样 时 出 现 缺失 。 为 了 克服 UV 覆盖 不 全 的 影响 ， 对 面 源 目标 的 空间 频率 信息 进行 更 好 的 采样 ， 可 以 通 
过 改变 子 卫星 和 主 卫星 的 相对 位 置 来 改变 系统 采样 频率 的 覆盖 范围 ， 获 得 之 前 缺失 的 空间 频率 采样 信 
息 ， 即 获得 不 同 基线 布局 下 的 干涉 图 。 再 经 过 图 像 合成 和 逆 滤 波 处 理 ， 最 终 得 到 包含 目标 较 完 整 信息 
的 图 像 。 通 过 这 样 的 方法 ， 使 合成 的 图 像 分 辩 率 等 效 为 一 个 更 大 口径 光学 系统 的 分 辨 率 ， 达 到 高 分 辨 
率 成 像 的 目的 。 

通常 的 合成 方法 是 将 这 些 在 不 同 基线 下 得 到 的 图 像 在 对 齐 的 基础 上 中 进行 琶 加 再 解 卷 积 ， 但 是 
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由 于 采集 的 图 像 受 子 望 远 镜 口径 较 小 、 积 分 时 间 受 限 等 各 种 因素 的 影响 ， 信 品 比 通常 都 不 高 。 仿 真实 
验 发 现 ， 当 得 到 的 干涉 图 在 信 噪 比较 低 的 情况 下 ， 采 用 直接 至 加 再 解 卷 积 的 处 理 方法 效果 不 理想 。 为 
了 提高 图 像 质量 ， 本 文 提 出 将 不 同 基 线 下 获得 的 图 像 通过 相应 的 频 域 滤波 器 ， 将 高 信 噪 比 的 频 域 信息 
提取 出 来 再 合成 ， 以 减少 合成 图 像 中 包含 的 噪声 ， 从 而 获得 更 好 的 合成 图 像 。 


1 原理 概述 


成 像 过 程 可 以 看 作 是 目标 物体 频率 函数 F(/.，/,) 与 光学 传递 函数 OTF Cf. f, ) 乘积 后 输出 像 频 谱 
函数 C(/.,/.) 。 在 非 相干 光照 射 下 ， 衍 射 受 限 系统 可 认为 是 只 改变 了 输入 的 各 频率 光 强 分 量 而 不 改 
变 它们 的 位 相 ， 光 学 传递 函数 在 数值 上 等 于 调制 传递 函数 MT; 
POF, Aff) ® POF, Mf) 

J[ Pes, 0d 


MTF(f,, f.) = |OTF(,, f) |= (1) 
其 中 ，@ 表 示 自 相关 符号 ; POf Mr 广 ) 是 孔径 函数 。 

稀疏 孔径 成 像 系统 在 理想 情况 下 可 以 看 作 是 由 N 个 直径 为 d 的 圆 形 孔径 构成 ， 在 数学 上 可 以 表 
示 为 一 个 子 孔径 光 瞳 函数 与 6 函数 阵列 的 卷 积 ， 即 


2 2 N 
P, eie [EET s Est ons (2) 
nz-zl 


Hub, PG, 站 ) 是 光 瞳 函数 ，N 是 子 孔径 个 数 ;cire( ) 是 圆 域 孔径 函数 ;(x ，y ) 是 第 ”个 子 孔 径 光 瞳 
的 中 心 坐标 。 

Golay3 是 结构 最 简单 的 二 维 阵列 结构 ， 由 3 个 圆 形 子 孔径 组 成 ， 而 且 这 3 个 圆心 依次 位 于 一 个 等 
边 三 角形 的 3 个 顶点 上 中。 

从 (2) 式 可 以 推导 出 Golay3 型 稀 朴 孔径 的 归 一 化 调制 传递 函数 为 


1 X; X3 
MTF = MTF, + —MTF, * ; .6 - l 一 3 
d 3 d DL, 6 Af , T] Af ). ( ) 


HB, (xx, yy) 表示 子 孔径 的 相对 位 置 。 

根据 (3) 式 ， 可 以 将 稀疏 孔径 系统 的 调制 传递 函数 看 作 是 由 分 布 在 不 同位 置 上 ， 强 度 为 1/3 倍 
( 零 频 除 外 ) 的 单个 子 孔径 的 调制 传递 函数 组 合 而 成 。 除 零 频 外 ， 其 它 调制 传递 函数 分 布 的 位 置 都 与 
子 孔径 排 布 有 关 ， 也 就 是 说 如 果 已 知 调制 传递 函数 的 位 置 分 布 就 可 以 推算 出 子 孔径 的 位 置 排 布 。 当 系 
统 调制 传递 函数 对 空间 采样 频率 上 覆盖 的 一 些 区 域 出 现 零 值 或 严重 衰减 时 ， 就 可 以 找到 相应 的 空间 位 置 
增加 和 孔径 排 布 进行 优化 ， 最 后 将 这 些 在 不 同 排 布 _ 
( 即 不 同 基线 ) 下 得 到 的 退化 图 进行 合成 处 理 。 A 

Jü.d LAS A de d ESI RO ET RS I. CCD pem io de P 
噪声 等 ， 总 噪声 分 布 近似 为 高 斯 分 布 ， 一 般 认 为 
是 零 均值 的 高 斯 白 噪 声 :5 ， 即 功率 谱 密 度 在 整 
个 频 域 内 均匀 分 布 的 噪声 。 已 知 单 幅 退化 图 的 噪 
声 功率 为 由， 如 果 将 退化 图 直接 丢 加 就 相当 于 直 
EAMT Y n, 倍 的 噪声 。 

通过 子 孔径 的 排 布 和 系统 调制 传递 函数 之 间 
的 计算 关系 ， 可 以 很 容易 了 解 到 退化 图 在 不 同 频 
域 范围 内 的 信 噪 比 情况 。 以 某 一 种 子 孔 径 排 布 得 
到 的 调制 传递 函数 为 例 ， 如 图 1。 图 中 ,中间 主 图 1 茶 子 孔径 排 布 的 调制 传递 函数 示意 图 
峰 的 峰值 表示 能 通过 系统 的 信号 强度 最 大 为 1 倍 ， Fig.1 The MTF of a sub-aperture configuration 
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周围 次 峰 则 表示 能 通过 的 信号 强度 最 大 只 有 原来 的 1/3， 其 余 值 为 0 的 部 分 ， 表示 相 应 频 域 范 围 内 无 
P P 
BARSE AEDA FA EA U HER IE LZ D 10182) , 1018 F) 0. 
Heb, POR LP. 分 别 代表 信号 和 噪声 的 有 效 功率 。 
因此 ， 本 文 针 对 各 种 基线 排 布 定义 了 相应 的 频 域 滤波 器 有 将 每 幅 图 中 频 域 信 噪 比较 高 部 分 的 
半 号 提取 出 来 再 合成 。 由 于 每 幅 干 涉 图 的 滤波 器 太 加 权 为 1， 而 且 范 围 不 完全 重复 ， 这 样 得 到 的 最 
终 合成 图 的 噪声 量 于 ,H,. n, 明显 小 于 通过 直接 希 加 得 到 的 ， 信 品 比 有 了 显著 的 提升 。 


2 ”图像 频率 提取 及 合成 方法 


设 单个 子 孔 径直 径 为 1m。 从 系统 子 孔 径 相 切 开始 ， 以 子 孔 径直 径 为 间距 向 外 挪动 子 孔 径 两 次 共 
得 到 3 种 不 同 的 排 布 。 此 时 外 接 最 大 包围 圆 直径 为 6. 15 m。 根 据 这 3 种 排 布 得 到 了 3 种 不 同 分 布 的 调 
制 传递 函数 。 由 于 发 现 空间 频率 信息 分 别 在 排 布 1 与 排 布 2、 排 布 2 与 排 布 3 截止 频率 交接 处 附近 明 
显 衰减 ， 所 以 增加 新 的 子 孔径 排 布 4 和 5， 对 这 两 交界 附近 的 调制 传递 函数 进行 补偿 和 增强 。 优 化 
后 ， 这 5 种 排 布 的 调制 传递 函数 在 最 大 截止 频率 方向 的 部 分 截面 示意 图 (在 同一 坐标 系 ) 如 图 2。 

与 稀 玻 孔径 系统 成 像 质 量 相当 的 单 孔 径直 径 称 作 等 效 口径 。 由 于 稀 下 孔径 系统 的 调制 传递 机 数 与 
单 孔 径 系统 相 比 ， 分 布 复 杂 ， 形 状 不 规则 ， 目 前 


! 尚 没有 计算 等 效 口径 的 统一 方法 四 。 从 图 2 可 以 
看 出 ， 当 前 系统 最 大 截止 频率 5.4/A(m) 明显 小 ; 
© 于 系统 包围 圆 所 对 应 的 6.15/Af(m) 。 所 以 选 最 。 os | 
一 大 截止 频率 对 应 的 口径 作为 系统 等 效 口 径 ， 此 处 F ost c Wis 
LÍ 为 5.4m。 2 
T 在 这 5 种 不 同 的 子 孔径 排 布 下 ， 获 取 5 幅 不 s 
- 同 的 图 像 ( 记 作 gi, as, ass g, Mgs) o HAX 5 幅 
图 像 通过 傅 里 叶 变换 到 频 域 ( 记 作 C, C, G, UA AAA 
G, 和 6 ) 。 根 据 各 自 调制 传递 函数 的 分 布 ， 分 (INE MEE UM MT RUE EE GE RECTO 
别 定 义 各 自 的 频 域 滤波 器 ， 如 (4) ~(8) 式 。 其 VAfim) 
e F, H, 为 低 通 滤 波 器 ; H,. H, H, VH, 分别 是 图 2 优化 后 各 排 布 的 调制 传递 函数 部 分 截面 图 
pr 不 同 通 带 范 围 的 带 通 滤波 器 。 Fig.2 The optimized sectional view of partial MTF 
6 l,rzR, for every sub-aperture configuration 
- H,(u, v) = (4) 
0,r>R 
1, R rR, 1, R Sr<R, 
H,(u, v)=10, r >R, , (5) H,(u, v) 10, r > R, ; (6) 
O,r«R, 0, r«R, 
1, R rR, 1, R Sr <R,;, 
H,(u, v) 2-30, r > R, ; (7) H,(u, v) 2-30, r > R; ; (8) 
0, r < R, 0, r< R, 


利用 这 些 频 域 滤波 器 可 以 将 5 幅 图 中 信 品 比较 好 的 频率 信息 提出 来 。 提 取 及 合成 方法 如 图 3。 图 
像 频 谱 合 成 方法 如 (9) 式 : 
Gan 7 CH, + GH, + GH, + GH, + G,H,. (9) 


频谱 直接 车 加 如 (10) 式 : 
GE 7 Cl * Ca * G, +G, * G,. (10) 
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-10 0 10 -10 0 10 -10 0 10 
f, Vf (m) f, Vf (m) 


(g) 


f, Mif (m) 


-10 0 10 
f, Vf (m) 


-10 0 10 -10 0 10 
f, Vf (m) f, Vf (m) f, Vf (m) 


图 3 空间 频率 信息 提取 及 合成 . (a) 原始 图 像 频率 信息 ; (b) 从 孔径 排 布 1 提取 的 频率 信息 ; (e). (d). (e) 
和 (f{) 分 别 对 应 从 孔径 排 布 2、3、4 de 5S 提取 的 频率 信息 ; (g) 合 成 图 像 频 率 信息 
Fig.3 Frequency extraction and synthesis. (a) The frequency of the original objective image; (b) The partial frequency extracted 
from the degraded image of configuration 1; (c) (d) (e) and (f) The corresponding partial frequency extracted from 


the degraded image of new configuration 2, 3, 4 and 5; (g) The final synthesized degraded image 


在 仿真 实验 中 引入 加 性 高 斯 白 噪 声 ， 当 信和 噪 比 为 35 dB 和 25 dB 时 分 别 得 到 退化 图 如 图 4。 


信 品 比 为 35 dB 时 : 


-200 -100 0 100 200 -200 -100 0 100 200 -200 -100 0 100 200 
信 噪 比 降 为 25 dB 时 : 
(c) 


-200 -100 0 100 200 -200 -100 0 100 200 -200 -100 0 100 200 
图 4 不 同方 法 得 到 的 仿真 图 像 . (a) 等 效 口 径 图 像 ; (b) 直接 登 加 得 到 的 合成 图 像 ; (e) 频 域 提取 方法 得 到 的 合成 图 像 
Fig.4  Degraded images obtained by different methods. (a) The degraded image obtained by effective Diameter method ; 
(b) The synthesized degraded image obtained by overlapping frequency directly; (c) The synthesized degraded 


image obtained by frequency filtrations 
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从 图 4 可 以 看 出 ， 当 图 像 信 品 比 下 降 后 ， 直 接 炙 加 的 合成 图 像 像 质 逐 渐变 差 ， 此 时 采用 频 域 提取 

的 合成 图 像 像 质 占 优 。 

3 图 像 复原 ( 解 卷 积 ) 及 评价 


由 不 同 基线 图 像 合成 的 图 像 需 经 过 解 卷 积 才能 更 好 地 恢复 图 像 。 逆 滤波 是 解 卷 积 最 直接 的 方法 ， 
设 G(u, vw) 为 退化 图 像 的 傅 里 叶 谱 ，H(u，, 5) 为 点 扩散 函数 的 傅 里 叶 谱 ， 则 估计 的 原始 图 像 的 傅 里 叶 
谱 u,v) 为 


Plu, v) = 一 一 (11) 


维 纳 滤波 (最 小 均 方 误差 滤波 ) 是 对 逆 滤 波 方法 的 改进 。 维 纳 滤 波 方法 综合 考虑 了 退化 函数 和 噪 
声 的 统计 特性 ， 以 减少 噪声 的 影响 ， 如 (12) 式 : 
x H'(u, v) 

|H(u, v)^|* y I5, (u, u)78(u, v)] ' 


Heb, RRR; S, (u, v)/S(Qu, o) 是 噪声 功率 谱 和 原始 信号 功率 谱 的 比值 。 当 噪声 的 统计 特性 未 知 时 ， 
常用 一 个 常数 天 近似 。 若 y=1，(w, o) 是 维 纳 滤波 器 ; 若 yY=0， 则 系统 变 成 了 单纯 的 去 卷 积 滤波 器 。 

由 于 图 像 经 过 频率 提取 和 合成 后 ， 合 成 图 的 品 信 比 有 较 明 显 的 起 伏 规律 。 因 此 ， 在 去 解 卷 积 滤波 
时 ， 对 维 纳 滤波 器 也 进行 了 相应 的 改进 : 在 对 应 不 同 的 区 域 ， 根 据 合 成 图 品 信 上 比 的 不 同 采用 了 不 同 的 
天 值 。 去 卷 积 后 的 复原 图 如 图 5, 


W(u, v) 


(12) 


信 噪 比 为 35 dB HJ : 
(a) 


-200 -100 0 100 200 -200 -100 0 100 200 -200 -100 0 100 200 


信 品 比 降 为 25 dB 时 : 


-200 -100 0 100 200 -200 -100 0 100 200 -200 -100 0 100 200 


图 5 不 同方 法 得 到 的 复原 效果 . (a) 等 效 满口 径 复原 图 ;(b) 直接 登 加 得 到 的 复原 图 ; 〈c) 频 域 提取 方法 得 到 的 合成 复原 图 


Fig. 5 ”Reconstructional images obtained by different methods. (a) The reconstructional image obtained by effective 


Diameter method; (b) The reconstructional synthetic image obtained by overlapping frequency directly ; 


(c) The reconstructional synthetic image obtained by frequency filtrations 
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对 于 图 像 的 复原 效果 ， 除 了 主观 判断 外 ， 采 用 传统 的 均 方 误差 ( Mean Squared Error, MSE) 法 和 当 
前 国际 上 比较 流行 的 SSIM (Structural Similarity) 法 。 均 方差 法 基于 误差 理论 ， 将 图 像 看 作 是 一 个 个 孤 
立 的 点 ， 通 过 计算 每 个 点 的 灰 度 值 差 异 而 对 图 像 进行 评测 ， 见 (13) 式 。SSIM 法 是 对 图 像 在 亮度 、 对 
比 度 、 结 构 3 个 层面 上 分 别 进行 量化 ， 通 过 加 权 相 乘 获得 结构 相似 值 ， 综 合 考虑 了 人 类 视觉 对 图 像 结 
构 性 信息 的 提取 功能 ， 比 均 方差 法 更 符合 人 类 的 视觉 特性 ， 见 (14) 式 : 


X 4Y 


2 5 1 ^ 2 
MSE = 4 = hm xy] [Us 0 -Aes y) Pidy, (13) 
MSSIMC, f) = SSIM = EIE, PG, forsan, PY. (14) 


其 中 , f(x, y) RR BOE IMS fx, y) 表 示 待 评价 图 像 ， a 、B、y 为 加 权 系 数 ; IF, f) cG D. 
sCf, 让) 为 相应 的 亮度 、 对 比 度 、 结 构 比 较 函 数 。 


均 方 根 值 越 小 ， 表 示 待 评价 图 像 和 参考 图 像 表 1 各 退化 图 客观 评价 计算 结果 
越 相 似 ; 结构 相 似 度 均值 越 大 ， 则 表示 结构 相 Table 1 Objective evaluation results for 
似 值 越 高 ， 说 明 待 评价 图 像 和 参考 图 像 的 相似 度 EE 
越 高 。 信 品 比 35 图 (a) 图 (b) BCe) 
在 不 同 信 噪 比 下 得 到 的 退化 图 经 过 解 卷 积 MSSIM 0.976 8 0. 7919 0. 800 3 
复原 后 ， 分 别 计算 均 方 根 和 结构 相似 度 均 值 ， 见 MSE 0. 000 6 0. 003 3 0. 003 1 
表 1。 信 噪 比 25 图 (a) 图 (b) Ele) 


um! 结合 图 5 和 表 1 可 以 看 出 ， 当 信 噪 比 为 35 dB 


ép CD MSSIM 0.916 6 0. 663 0 0.744 1 
im i dn - oo gA MSE 0.0019 0. 003 5 0. 003 0 

(b) 的 客观 评价 指标 只 是 稍微 好 一 点 ， 从 主 和 

图 (a) 是 以 原 图 为 评价 参考 , 图 (b) 和 图 (c) 都 是 以 图 (a 

观 上 观察 ， 人 的 眼睛 基本 分 不 出 哪 幅 图 更 好 。 2 Mida 


然而 当 信 品 比 下 降 到 25 dB 时 ， 对 合成 图 (e) Fg (a) is with reference to the original image; (b) and (c) 
客观 评价 指标 就 明显 比 对 直接 琶 加 图 像 (D) 的 客 ue esce to e (i 
观 评 价值 要 好 一 些 。 在 用 人 眼 进行 主观 评价 时 ， 
能 明显 看 出 合成 图 (ec) 中 的 图 像 细 节 更 清晰 。 


O 4 结论 


当 原始 图 像 的 信 噪 比 受 限时 ， 利 用 频 域 抽取 算法 将 在 不 同 子 孔 径 排 布 下 得 到 的 干涉 图 合成 为 一 幅 
包含 较 完 整 目标 信息 的 退化 图 像 ， 再 经 过 改进 的 维 纳 滤波 算法 得 到 了 最 终 较 清晰 的 复原 图 像 。 由 于 保 
留 了 各 幅 干 涉 图 频 域 中 信 噪 比较 高 的 部 分 ， 在 解 卷 积 复原 后 ， 其 图 像 质量 比 直接 县 加 得 到 的 有 明显 提 
升 ， 采 用 人 眼 判 读 和 客观 评价 的 方法 都 显示 ， 该 方法 能 更 好 地 恢复 原 图 像 的 结构 信息 。 
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Abstract: A Fizeau-typesparse-aperture imaging system made of a constellation of telescope-bearing 
satellites can potentially achieve extremely high resolution. But in such a system, the Modulation Transfer 
Function ( MTF) usually contains a zero value. This means that the system fails to acquire a particular spatial 
frequency. To overcome this problem, we acquire multiple degraded images by repeatedly changing the sub- 
aperture configuration and then synthesize their spatial frequencies to create a single improved image. We first 
analyze the MTFs between each individual sub-aperture and the whole sparse aperture system. Then, we design 
new sub-aperture configurations that acquire more spatial frequency and use frequency-domain filters to extract 
the high-SNR part from each degraded image. Finally, we synthesize these frequencies before image 
reconstruction. Our results demonstrate that this method effectively improves the quality of the final synthesized 
image by obtaining as much spatial frequency information from all degraded images as possible. 


Key words: Sparse aperture; Image reconstruction; Frequency-domain filtering; Image synthesis 


